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I.  INTRODUCTION 


This  paper  is  a  sequel  to  an  earlier  work  on  muzzle  flash  and  its 
suppression  from  an  elementary  kinetic  point  of  view.  In  the  earlier  report 
we  discussed  definitions  of  terras,  the  importance  of  the  problem,  the  options 
available  to  suppress  the  flash,  and  several  apparently  contradictory 
laboratory  results  on  flame  inhibition.  The  present  paper  builds  upon  the 
recent  work  of  Hynes,  Steinberg,  and  Schofield2  to  produce  a  hypothetical  but 
realistic  suppression  network  of  elementary  chemical  reactions.  The  most 
important  aspect  of  the  present  work  is  the  methodology  employed  to  construct 
the  network.  The  methodology  is  quite  general  and  can  be  used  iteratively  as 
more  accurate  data  (species  or  rate  coefficients)  become  known. 


We  shall  begin  with  a  brief  review  of  the  problem  (see  Reference  1  for 
more  details),  then  illustrate  the  methodology  by  actually  constructing  a 
suppression  network. 


The  effluent  from  a  large  artillery  piece  is  known  to  consist  of  a  fuel 
rich  mixture  containing  H2  and  CO.  Upon  mixing  with  the  surrounding  air,  a 
combustion  process  can  take  place  that  results  in  a  high  intensity  flash. 
This  phenomenon  is  called  secondary  muzzle  flash  or  simply  muzzle  flash. 
About  one  half  the  energy  of  the  propellant  can  to  be  found  in  this  effluent, 
so  that  muzzle  flash  can  be  accompanied  by  a  substantial  blast.  In  the  past 
both  this  flash  and  blast  have  been  suppressed  in  many  weapons  by  the  addition 
of  alkali  salts,  usually  potassium  sulfate  or  potassium  nitrate.  Since  only  a 
few  percent  (by  weight)  salt  addition  has  been  found  to  be  effective  in  many 
gun  systems,  the  detailed  mechanism  appears  to  be  chemical  rather  than 
physical  in  nature.  If  this  mechanism  were  understood,  then  measures  to 
suppress  the  flash  could  be  taken  into  account  early  in  the  design  stages  of  a 
new  weapon  or  weapon  subsystem.  This  procedure  would  lead  to  substantial  cost 
saving  over  the  current  practice  of  treating  each  case  empirically  after  the 
fact. 


Earlier*  we  discussed  the  problems  that  surround  a  popular  suppression 
network  that  consists  of  two  reactions.  Recent  work2  has  shown  how  complex 
such  a  network  can  be.  In  this  paper  we  shall  build  upon  this  work  and,  in  a 
systematic  fashion,  construct  an  advanced  suppression  network  that  forms  a 
working  hypothesis  against  which  future  experiments,  both  in  the  field  and  in 
the  laboratory,  can  be  tested. 


1.  J.  M.  Heimerl,  "Muzzle  Flash  Suppression  and  Alkali  Salt  Inhibition," 
JANNAF  13th  Plume  Technology  Meeting,  Houston,  TX,  Vol  I,  pp.  413-418,  27-29 
April  1982.  Also  available  as:  J.  M,  Heimerl,  "Muzzle  Flash  and  Alkali  Salt 
Innibition  from  an  Elementary  Kinetic  Point  of  View,"  Technical  Report  ARBRL- 
TR-02479,  March,  1983.  AD  A126  129. 

2.  A.  J.  Hynes,  M.  Steinberg  and  K.  Schofield,  "The  Chemical  Kinetics  and 
Thermodynamics  of  Sodium  Species  in  Oxygen-Rich  Hydrogen  Flames,"  J^.  Chem. 
Ph v s . .  Vol. 80,  pp.  2  385-2597  ,  1  984. 
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II.  THE  CONSTRUCTION  OF  THE  SUPPRESSION  NETWORK 

As  with  any  new  reaction  network,  three  things  must  be  determined:  the 
reaction  species,  the  network  of  reactions(  and  the  rate  coefficients  for  these 
reactions. 

A.  THE  SPECIES 

The  species  that  correspond  to  hydrogen  as  a  fuel  are:  H,  0,  OH,  Hj,  O2, 
H02,  and  H2O.  To  include  carbon  monoxide  as  a  fuelfthe  species  CO  and  CO2 
must  be  added  to  this  list.  Since  we  are  dealing  with  air  we  need  also 
consider  N2.  However,  the  temperatures  are  such  that  *2  act8  only  as  a 
diluent.  The  species  H2O2  is  absent  from  the  above  list.  The  reason  for  this 
neglect  is  as  follows.  The  major  formation  path  for  H2O2  is  given  by  the 
three-body  reaction: 

OH  ♦  OH  ♦  M  ■  H202  +  M, 

whose  reaction  rate  coefficient  is  given^  by:  1.3  x  10^2  x  T“2  cm®/mole2-sec. 
At  one  atmosphere  and  a  temperature  of  1500  K  the  effective  two-body  rate 

coefficient  for  this  reaction  is  about  5  x  10*®  cm^/mole-sec.  Two-body 
reactions  involving  the  OH  radical  have  much  greater  values  for  their  rate 
coefficients  at  1500  K.  For  example,  the  reactions: 

OH  ♦  0  »  02  ♦  H, 

OH  ♦  H2  *  H2O  +  H,  And 

OH  ♦  OH  -  H20  ♦  0 

have  values  of  about  2  x  10*^,  1  x  10*2  and  6  x  10*2  cm^/raole-sec, 
respectively.-*  Thus,  the  species  H2O2  is  not  considered  because  at  high 
temperatures  its  formation  rate  is  much  too  small  compared  to  other  reactions 
involving  OH. 

There  remains  the  selection  of  the  suppressant  species.  The  species  Na, 
NaO,  NaO^  and  NaOH  were  required2  to  describe  adequately  the  inhibition  of 
oxygen-ricK  hydrogen  flames  by  sodium  salt.  By  analogy  we  consider  the 
corresponding  potassium  molecules:  K,  KO,  K02  and  KOH.  In  addition,  since  the 
muzzle  effluent  is  fuel-rich  we  also  consider  the  species  KH. 

Thus  we  shall  consider  14  reactive  species  and  the  diluent  N2> 


3.  J.  Warnatz,  "Survey  of  Rate  Coefficients  in  the  C/H/0  System,"  Sandia 
Report  SAND83-8606,  February  1983. 
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Tablt  shows  r'  '  i.;t  of  14  active  species,  their  standard  enthalpies  of 
f  orma  l  ion,  and  tht  ir  absolute  entropies  (both  at  298  K).  These  values  are 
taken  froi.i  standard  sources,  except  tor  the  enthalpy  of  KO2,  which  is 
estimated  in  the  appendix.  Conversion  to  SI  units  can  be  made  through  the 
relation  4.1b4  joules  -  1  calorie. 

TABLE  1.  LIST  OF  ACTIVE  SPECIES 


enthalpy  entropy 

(kcal/mole)  ( eal/mole/K) 


1 

H 

52.10 

27.40 

a  (p. 

106) 

2 

K 

21.30 

38.30 

b 

3 

0 

59.60 

38.50 

a  (p. 

132) 

4 

Ha 

0.00 

31.20 

b 

5 

KH 

26 .30 

47.30 

c,  b 

6 

K0 

17.50 

56.90 

c,  b 

7 

OH 

9.50 

43.90 

d 

8 

°2 

0.00 

49.00 

b 

9 

ho2 

3.50 

54.40 

e,  b 

10 

HoO 

-57.80 

45.10 

b 

11 

KOH 

-54.50 

56.90 

c,  b 

12 

K07 

-25.00 

62.00 

f 

13 

CO 

-26.40 

47.20 

b 

14 

co2 

-94.10 

51.10 

b 

a ) 

J.  Phvs.  &  Chem.  Ref. 

Data.  Vol. 4.  1975. 

b> 

JANAF 

Thermochemical 

Tables,  Second  Edition, 

1971. 

c) 

J.  Chem.  Thermo..  Vol 

.14,  pp. 1103-1113,  1982 

• 

d) 

Ls  Phy 

rs.  6  Chem.  Ref. 

Data.  Vo  1.3,  p.443,  1974. 

e) 

Is. 

r s.  Chem..  Vol. 87,  pp. 3479-3482 ,  1983. 

f) 

This  work,  see  append 

ix. 

B.  IHE  KEACT10N  NETWORK 

The  next  step  in  this  process  is  the  construction  of  a  complete  reaction 
network.  We  consider  reactions  involving  as  many  as  three  molecules;  i.e.,  we 
neglect  reactions  with  four  or  more  collision  partners.  We  want  to  be  sure 
that  we  have  considered  them  all  since  the  neglect  of  one  implicitly  assumes 
that  its  rate  coefficient  is  zero.  In  principle,  one  simply  writes  down  all 
combinations  of  the  species  to  form  the  reactions.  However,  doing  this 
manually  for  more  than  two  or  three  species  proves  to  be  tedious  and  subject 
to  error.  Therefore,  following  Indritz,  Maday  and  Sheinson,  vs  have 
constructed  a  computer  program  to  perform  the  necessary  bookkeeping  operations 
to  assure  that,  nil  two-  and  three-body  reactions  are  in  fact  considered. 
Since  the  forward  and  reverse  directions  of  an  elementary  reaction  are 


4.  D.  Indritz,  M.  Maday  and  R.  Sheinson,  "Construction  of  Large  Reaction 
Mechani  .,n,s,"  NHL  Report  8498,  October  1981. 
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automatically  taken  into  account  and  since  the  reverse  of  three-body  reactions 
of  the  form  A  ♦  B  ♦  M  <->  AB  +  M  are  unioolecuiar  reactions,  we 
automatically  obtain  the  unimolecular  decomposition  reactions. 

Table  2  shows  the  list  of  all  exothermic  two-body  reactions  that  can  be 
made  from  the  species  listed  in  Table  1.  In  order  to  avoid  duplication,  only 
those  reactions  that  are  exothermic  at  298  K  (i.  e. ,  have  negative  Gibbs  free 
energies)  are  printed  out.  The  selection  of  the  Gibbs  energy  rather  than  the 
enthalpy  of  the  reaction  is  arbitrary. 

We  can  now  begin  the  systematic  elimination  of  the  less  important 
reactions.  We  make  the  reasonable  assumption  that  if  more  that  one  bond  is 
broken  in  a  given  reaction,  it  will  have  a  very  small  rate  coefficient  and  so 
can  be  neglected.  With  the  exception  of  KO2,  the  structures  of  all  the 
species  listed  in  Table  1  are  known.  We  assume  that  the  gas  phase  KO2 

molecule  has  the  same  triangular  geometry  as  the  solid  phase  CaC2  molecule.^ 
This  assumption  is  consistent  with  our  determination  of  the  heat  of  formation 
of  K(>2>  (See  appendix.) 

We  find  that  reactions  19,  22,  24-26,  30,  32-39,  47,  52  and  54  require 
more  than  one  bond  to  be  broken  and  so  are  deleted  from  the  network.  Note 
that  reactions  25,  26,  34,  38f  and  52  contain  KO2  but  require  more  than  one 
bond  to  be  broken  no  matter  what  the  KO2  structure.  On  the  other  hand, 
reactions  8,  10,  49, and  55  do  depend  on  the  structure  assumed  for  KO2.  With  a 
triangular  structure,  these  reactions  require  more  than  one  bond  to  be  broken 
and  so  they  too  are  eliminated  from  the  network. 

As  shown  in  Table  3,  all  reactions  can  conveniently  be  classified  into: 
those  reactions  that  are  involved  in  the  description  of  the  hydrogen-oxygen 
flame,  those  reactions  that  describe  the  conversion  of  carbon  monoxide  to 
carbon  dioxide,  those  reactions  that  are  the  potassium  analogs  of  the  sodium 
network  of  reference  2,  those  reactions  that  involve  the  species  KH,  and  those 
reactions  deleted  (two  categories). 

Table  4  lists  three-body  reactions  generated  by  the  program.1  Not  shown 
in  this  Table  are  the  reactions: 

c 

0  +  Hj  ♦  M  ■  H2O  +  M,  and 

0  +  KH  +  M  -  KOH  +  M,  ! 

c  4  c  ;  .  1  c.. 

because  it  is  thought  that  the  required  rearrangement  to  form  the  product 
molecule  is  improbable.  In  addition,  the  four  reactions: 


5.  For  a  discussion  of  the  structure  of  these  molecules  in  the  solid  phase, 
see  for  example,  ADVANCED  INORGANIC  CHEMISTRY.  A  Comprehensive  Text,  by  F.  A. 
Cotton  and  G.  Wilkinson,  Wiley  and  Son,  New  York,  fourth  edition,  p.  363  and 
p.  498,  1980. 


. ^RLE  >.  TWO-BODY  REACTIONS 


reac  t ion 


1 

H 

t  KH 

=  K 

4  ll2 

2 

H 

+  KO 

=  K 

+  OH 

3 

H 

+  OH 

»  0 

4-  H2 

4 

H 

+  HO, 

»  0 

+  ll20 

5 

H 

+  HO-j 

=  h2 

4-  02 

6 

H 

+  H(;“ 

=  OH 

4  OH 

7 

H 

+  Kofi 

=  K 

4  H20 

8 

9 

H 

H 

+  Kn-, 

+  ko2 

=  K 
=  0 

4  H02 
4  KOH 

10 

H 

+  KO-) 

»  KH 

4  02 

11 

H 

+  KOj 

=  KO 

4  OH 

12 

K 

+  !:«  2 

=  0 

4  KOH 

13 

K 

+  H02 

=  KO 

4  OH 

14 

0 

+  KH 

-  H 

4  KO 

15 

0 

+  KH 

-  K 

4  OH 

16 

0 

+  KC 

=  K 

4  02 

17 

0 

+  OH 

-  H 

+  02 

18 

19 

20 

0 

0 

0 

+  iio2 
+  H;0 
♦  K00 

-  OH 

-  h2 
«  KO 

+  02 
+  Oj 

4  On 

21 

h2 

+  Kt" 

«  H 

♦  KOH 

22 

h2 

+  KO 

-  K 

4  H20 

23 

h2 

♦  OH 

-  H 

4  H20 

24 

h2 

+  HOj 

-  OH 

4  H20 

25 

h2 

+  ko2 

=  KO 

4  H20 

26 

h2 

+  KO, 

«  OH 

4  KOH 

27 

KH 

+  KC*" 

=  K 

4  KOH 

28 

KH 

+  OK 

**  H 

4  KOH 

29 

KH 

+  Oil 

*  K 

4  H20 

30 

KH 

+  OH 

-  h2 

4  KO 

31 

KH 

+  0  , 

-  K 

4  H02 

32 

KH 

+  0* 

=  0 

4  KOH 

33 

KH 

+  On 

=  KO 

4  OH 

34 

KH 

+  h62 

•  h2 

4  K02 

35 

KH 

+  ho2 

=  KO 

4  H20 

36 

KH 

+  Ho  2 

=  OH 

4  KOH 

37 

KH 

+  h2o 

=  h2 

4  KOH 

38 

KH 

+  K02 

=  KO 

4  KOH 

39 

KH 

*■  co2 

=  KOH 

4  CO 

40 

KO 

+  KO 

=  K 

4  K02 

41 

KO 

+  OH 

=  0 

>  KOH 

42 

KO 

+  H07 

=  OH 

4  K02 

43 

KO 

+  mb 

=  o2 

4  KOH 

44 

KO 

H-,0 

=  OH 

4  KOH 

45 

KO 

+  CO 

=  K 

4  C02 

46 

OH 

+  OH 

=  0 

+  H20 

47 

OH 

+  OH 

=  T ' 

2 

4  02 

ent  ha Idv 

entroDY 

free  enerev 

-57.1000 

-5.2000 

-55.5504 

1 

-38.8000 

-2.1000 

-38.17A2 

2 

-2.0000 

-1 .6000 

-1.5232 

3 

-53.8000 

1.8000 

-54.3364 

4 

-55.6000 

-1.6000 

-55.1232 

5 

-36.6000 

6.0000 

-38.3880 

6 

-34.1000 

-0 . 9000 

-33.8318 

7 

-2.3000 

3.3000 

-3.2834 

8 

-22.0000 

6.0000 

-23.7880 

9 

-0.8000 

6.9000 

-2.8562 

10 

-0.1000 

11.4000 

-3.4972 

11 

-19.7000 

2.7000 

-20.5046 

12 

2.2000 

8.1000 

-0.2138 

13 

-16.3000 

-1.5000 

-15.8530 

14 

-55.1000 

-3.6000 

-54.0272 

15 

-55.8000 

-8.1000 

-53.3862 

16 

-17.0000 

-6.0000 

-15.2120 

17 

-53.6000 

0.0000 

-53.6000 

18 

-1 .8000 

-3.4000 

-0.7868 

19 

-17.1000 

5.4000 

-18.7092 

20 

-19.9000 

-3.8000 

-18.7676 

21 

-54.0000 

-4.7000 

-52.5994 

22 

-15.2000 

-2.6000 

-14.4252 

23 

-51.8000 

3.4000 

-52.8132 

24 

-15.3000 

8.8000 

-17.9224 

25 

-20.0000 

7.6000 

-22.T648 

26 

-77.0000 

-9.0000 

-74.3180 

27 

-38.2000 

-6.9000 

-36.1438 

28 

-72.3000 

-7.8000 

-69.9756 

29 

-18.3000 

-3.1000 

-17.3762 

30 

-1.5000 

-3.6000 

-0.4272 

31 

-21 .2000 

-0.9000 

-20.9318 

32 

0  7000 

4.5000 

-0.6410 

33 

-54.8000 

-8.5000 

-52.2670 

34 

-70.1000 

0.3000 

-70.1894 

35 

-74.8000 

-0.9000 

-74.5318 

36 

-23.0000 

-4.3000 

-21.7186 

37 

-38.3000 

4.5000 

-39.6410 

38 

-13.1000 

5 . 7000 

-14.7986 

39 

-38.7000 

-13.5000 

-34.6770 

40 

-21.9000 

-5.4000 

-20.2908 

41 

-36.5000 

-5.4000 

-34.8908 

42 

-75.5000 

-5.4000 

-73.8908 

43 

-4.7000 

-1.2000 

-4.3424 

44 

-63.9000 

-14.7000 

-59.5194 

45 

-17.2000 

-4.2000 

-15.9484 

46 

-19.0000 

-7.6000 

-16.7352 

47 
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TABLE  2 


TWO-BODY  REACTIONS  (continued) 


reaction 


48 

OH 

+ 

HO, 

3? 

°2 

h2o 

49 

OH 

+ 

KO, 

a 

°2 

+ 

K0H 

50 

OH 

CO 

m 

H 

+ 

co2 

SI 

°2 

+ 

CO 

9 

0 

♦ 

co2 

52 

ho2 

+ 

K0H 

9 

H,0 

+ 

ko2 

53 

ho2 

+ 

CO 

a 

OH 

•f 

co2 

54 

h2o 

+ 

CO 

- 

H, 

♦ 

co2 

55 

ko2 

♦ 

CO 

a 

ICO 

+ 

co2 

enthalDV 

entroDV 

free  enerzv 

-70.8000 

-4.2000 

-69.5484 

48 

-39.0000 

0.0000 

-39.0000 

49 

-25.1000 

-12.6000 

-21.3452 

50 

-8.1000 

-6.6000 

-6.1332 

51 

-31.8000 

-4.2000 

-30.5484 

52 

-61.7000 

-6.6000 

-59.7332 

53 

-9.9000 

-10.0000 

-6.9200 

54 

-25.2000 

-1.2000 

-24.8424 

55 

TABLE  3.  CATEGORIES  OF  REACTIONS 


classification 


reaction  numbers 


h2/o2 

co/co2 

K  analogs  to  Na  scheme 
reactions  involving  KH 


3-6,17-18,23,46,48 

45,50-51,53 

2 , 7 , 9 , 1 1-13 , 16 , 20-21 , 40-44 
1,14-1!, 27-29,31 


deleted  -  multiple  bond 

breaking  19,22 ,24-26 ,30,32-39,47 ,52,54 

deleted  -  depends  on 

K02  structure  8,10,49,55 


TABLE  4. 

THREE- BODY 

REACTIONS 

H 

c  c- 

f»2 

<. 

K 

KH 

H 

+  ' 

KO 

♦ 

M 

KOH  + 

M 

°2 

ho2 

* 

O 

KO 

K 

OH 

+ 

M 

< 

KOH  + 

M 

°2 

ko2 

J 

(ko 

i 

fKOz 

O 

*  1 

[co 

+ 

M 

I 

^co2  + 

M 

12 
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h  +  0  t  !!  -  OH  *  K , 

II  •  OH  v  M  =  H 2 0  +  M, 
o  «•  0  -*•  ll  -  (>2  +  M,  and 

0  +  OH  *  M  =  I102  ♦  M, 

do  not  appear  because  competing  two-body  reactions  involving  these  radicals 
are  known  to  be  more  important. 

C.  THE  RATE  COEFFICIENTS 


There  are  several  compilations  and  evaluations  of  rate  coefficients  for 


the  hydrogen-oxygen  network^ 
recommendations  of  Warnatz.' 


3,6,7 


and  ve  have  elected  to  use  the  most  recent 


He  point.,  out  that  reaction  4  of  Table  2: 

H  +  H02  -  H20  +  0, 

has  been  postulated,  even  though  there  is  no  direct  evidence  of  its 
occurrence.  In  addition,  he  argues  that  indirect  determinations  show  its  rate 
coefficient  to  be  small  relative  to  the  other  two  branches  of  the  H  +  H02 
reaction  (reactions  5  and  6  of  Table  2).  We  accept  his  arguments  and  so  do 
not  consider  this  reaction  further. 

Warnatz-*  has  also  reviewed  the  values  for  the  C0/C02  rate  coefficients 
and  we  have  again  used  his  recommendations.  The  one  reaction  of  the  C0/C02 
group  not  reviewed  in  the  literature  is: 

K0  +  CO  =  K  +  C02. 

Fortunately,  it  is  known  that  the  sodium  analog  of  this  reaction  is  fast  and 

so  we  assume  a  value  of  1  x  10^  cm^/mole-sec  for  its  rate  coefficient.  We 
use  the  published  values  of  the  rate  coefficients  for  the  sodium  network  as 
those  values  appropriate  to  the  potassium  analog  reactions.  (See  Table  3.) 
The  single  exception  to  this  strategy  is  the  value  of  the  rate  coefficient  for 
the  reaction: 


K  +  02  ♦  M  =  K02  +  M. 


6.  G.  Dixon-1  ewis,  "Kinetic  Mechanism,  Structure  and  Properties  of  Premixed 
Flames  in  Ily i‘ r og e n-Oxy g en-N i t r o gen  Mixtures,"  Proc.  Rov.  Soc.  London. 
Vo  1.2  92  ,  pp.4  3-99,  1  97  9. 

7.  N.  Cohen  and  K.  R.  Westberg,  "Chemical  Kinetic  Data  Sheets  for  High 
Temperature  Chemical  Reactions."  J.  Phvs.  Chem.  Ref.  Data.  Vol.12,  pp.  531— 
566  ,  1  983. 

8.  J.  Silv.i,  Aerodyne  Research  Inc.,  private  communication,  1984. 
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Here  we  have  adopted  Husain  and  Plane's^  value  for  the  corresponding  sodium 

reaction,  3  x  lO1^  cm^/mole-sec,  independent  of  temperature. 

Those  two-body  rate  coefficients  for  which  we  could  find  no  published 

values  were  set  at  1  x  lO1^  cm^/mole-sec  (  ■  0.17  x  10-0^  cm^ /par tide- sec). 
Similarly  the  three-body  rate  coefficients  were  assumed  to  have  a  rate 

coefficients  that  varied  as  A  x  T-1 ,  where  A  *  lO1^  cm^/mole^-sec 

a  n  A 

(  ■  0.83  x  10“JU  cmJ/particle  -sec).  This  network  of  48  reactions  and  their 
respective  rate  coefficient  parameters  is  shown  in  Table  5. 


III.  RESULTS  AND  DISCUSSION 

At  this  point,  the  network  of  reactions  shown  in  Table  5  is  ready  for 
testing.  The  MEFF  muzzle-flash-prediction  code10  is  the  only  operational 
flash  prediction  code  with  detailed  chemistry,  it  has  been  compared  with  other 
available  codes  and,  to  some  extent,  .with  observations,11  and  a  user  guide  has 
been  written  to  facilitate  its  use.1*  Reference  11  documents  calculations 
made  with  MEFF  using  an  earlier,  less  comprehensive  25-reaction  network.  We 
found  that,  in  a  nominal  155-mm  howitzer  simulation  with  a  nominal  M203 
propelling  charge,  1Z  suppressant  permitted  flash,  but  2.7Z  K2SO4 
suppressant  suppressed  the  flash.  Similarly,  in  an  81-mm  mortar  simulation,  a 
nominal  M10  propellant  lacking  suppressant  flashed,  but  the  addition  of  2Z 
K2SO4  suppressed  the  flash,  in  agreement  with  observations. 

We  repeated  the  howitzer  and  mortar  calculations  with  the  48-reaction 
network  of  Table  5,  and  the  predictions  of  the  MEFF  flash-prediction  code  were 
nearly  identical. 

There  still  remained  the  question  of  whether  or  not  the  network  of  48 
reactions  was  a  minimum  needed  to  describe  flash  suppression.  With  the  goal 
now  of  finding  out  whether  the  size  of  this  reaction  network  could  be  further 
reduced,  we  took  reactions  1-41  of  Table  5  as  the  reaction  network  to  describe 
the  inhibition  of  a  hydrogen-air  flame.  We  employed  our  steady-state, 


9.  D.  Husain  and  J.  M.  C.  Plane,  "Kinetic  Investigation  of  the  Reaction 
between  Na  +  Oj  +  M  by  Time-Resolved  Atomic  Resonance  Absorption 
Spectroscopy."  J.  Chem.  Soc.  Faraday  Trans.  2.  Vol.78,  pp.  163-178,  1982. 

10.  V.  Yousefian,  "Muzzle  Flash  Onset,"  Contractor  Report  ARBRL-CR-00477 , 
February,  1982.  AD  B063  573L. 

11.  G.  E.  Keller,  "An  Evaluation  of  Muzzle  Flash  Prediction  Models,"  U.  S. 
Army  Ballistic  Research  Laboratory  ARBRL-MR-03318,  November  1983. 

12.  G.  E.  Keller,  "A  MEFF  User's  Guide,"  U.  S.  Army  Ballistic  Research 
Laboratory  ARB RL- MR-03 3 62,  July  1984. 
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TABL..  5.  NETWORK  OF  .V.L  SINGLE-BOND-BREAKING  REACTIONS 


reaction 

A 

B 

c 

1  OH 

+  »2 

-  h2o 

+  1( 

0.17E-15 

-1.6 

-3298.8 

2  H 

+  o2 

-  OH 

+  0 

0.20E-06 

0.9 

-16573.2 

3  0 

+  H2 

=  OH 

+  H 

0.25E-16 

-2.0 

-7551.4 

4  H 

+  02  +  M 

=  ho2 

+  M 

0.55E-29 

0.8 

0.0 

5  H 

♦  ho2 

-  OH 

+  OH 

0.25E-09 

0.0 

-1003.5 

6  H 

+  ho2 

“  H2 

+  o2 

0.42E-10 

0.0 

-693.5 

7  OH 

+  ho2 

=  h2o 

+  o2 

0.33E-10 

0.0 

0.0 

8  0 

+  ho2 

-  OH 

+  0, 

0.33E-10 

0.0 

0.0 

9  H 

+  H  +  M 

’  H2 

+  M 

0.50E-29 

1.0 

0.0 

10  H 

+  H  +  H2 

*  h2 

+  h2 

0.27E-30 

0.6 

0.0 

11  OH 

+  OH 

-  0 

♦  h2o 

0.25E-14 

-1.1 

0.0 

12  H 

♦  ko2 

=  KO 

♦  OH 

0.12E-09 

0.0 

0.0 

13  KO 

+  h2o 

-  OH 

+  KOH 

0.17E-09 

0.0 

0.0 

14  H 

+  KOH 

-  K 

+  h20 

0.17E-10 

0.0 

0.0 

15  H 

+  KO 

-  K 

+  OH 

0.33E-10 

0.0 

0.0 

16  H 

+  ko2 

-  0 

+  KOH 

0.33E-11 

0.0 

0.0 

17  H 

♦  K02 

-  K 

+  ho2 

0.33E-10 

0.0 

0.0 

18  K 

♦  ho2 

-  0 

♦  KOH 

0.66E-13 

0.0 

0.0 

19  K 

+  ho2 

-  KO 

+  OH 

0.33E-11 

0.0 

0.0 

20  0 

+  KO 

-  K 

♦  o2 

0.33E-11 

0.0 

0.0 

21  0 

+  ko2 

-  KO 

+  o2 

0.17E-11 

0.0 

0.0 

22  H2 

+  KO 

-  H 

♦  KOH 

0.83E-09 

0.0 

0.0 

23  KO 

+  OH 

■  0 

*  KOH 

0.10E-09 

0.0 

0.0 

24  KO 

♦  ho2 

-  OH 

+  ko2 

0.33E-11 

0.0 

0.0 

25  KO 

+  ho2 

-  o2 

♦  KOH 

0.17E-09 

0.0 

0.0 

26  OH 

+  ko2 

"  °2 

+  KOH 

0.50E-12 

0.0 

0.0 

27  K 

+  OH  +  M 

*  KOH 

+  M 

0.25E-30 

1.0 

0.0 

28  K 

+  02  +  M 

-  ko2 

+  M 

0.83E-30 

0.0 

0.0 

29  KO 

+  KO 

-  K 

+  KO, 

0.17E-09 

0.0 

0.0 

30  H 

♦  KO  +  M 

-  KOH 

+  M  1 

0.83E-30 

1.0 

0.0 

31  0 

+  KO  +  M 

-  KO, 

+  M 

0.83E-30 

1.0 

0.0 

32  K 

♦  0  +  M 

-  KO 

+  M 

0.83E-30 

1.0 

0.0 

33  H 

+■  KH 

-  K 

+  h2 

0.17E-09 

0.0 

0.0 

34  0 

♦  KH 

-  H 

+  KO 

0.83E-10 

0.0 

0.0 

35  0 

+  KH 

-  K 

+  OH 

0.83E-10 

0.0 

0.0 

36  KH 

♦  KO 

-  K 

+  KOH 

0.17E-09 

0.0 

0.0 

37  KH 

+  OH 

-  H 

+  KOH 

0.17E-10 

0.0 

0.0 

38  KH 

+  OH 

-  K 

♦  H20 

0.17E-09 

0.0 

0.0 

39  KH 

+  o2 

-  K 

♦  ho2 

0.17E-09 

0.0 

0.0 

40  H 

+  ko2 

-  KH 

♦  0, 

0.17E-09 

0.0 

0.0 

41  H 

+  K  +  M 

»  KH 

+  M 

U.83E-30 

1.0 

0.0 

42  CO 

+  OH 

-  H 

+  co2 

0.73E-17 

-1.5 

741.0 

43  CO 

+  o2 

=  0 

+  co2 

0.42E-11 

0.0 

-47800.0 

44  CO 

+  ho2 

«  OH 

+  co2 

0.25E-09 

0.0 

-23590.0 

45  CO, 

+  0 

=  CO 

+  o2 

0.28E-10 

0.0 

-52580.0 

46  C02 

+  H 

*  CO 

+  OH 

0.27E-09 

0.0 

-26300.0 

47  0 

+  CO  +  M 

=  CO, 

+  M 

0.15E-33 

0.0 

4540.0 

48  KO 

♦  CO 

»  K 

+  co2 

0.17E-09 

0.0 

0.0 

The  rate 

coefficient 

■  A  x  T 

-B  x  eC/RT 

where  A  in  units 

of  cm-particle- 

sec 

C  in  units  of  calorie/mole 
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laminar,  flat-flame  model13,14  as  a  tool  to  make  these  decisions.  The 
selection  of  input  parameters  other  than  the  reactions  and  their  rate 
coefficients  will  now  be  outlined. 

The  transport  and  thermodynamic  parameters  are  shown  in  Table  6.  The 
first  line  identifies  the  species  and  the  five  transport  parameters.  The 
first  is  the  mass  of  the  species,  the  second  and  third  are  the  Lennard-Jones 
collision  diameter  (in  Angstroms)  and  the  Lennard-Jones  potential  well  depth 
(in  Kelvins),  respectively,  the  fourth  rrnmber  is  the  dipole  moment  (in  Debyes, 
where  one  Debye  ■  1 0"1 ® c m 3 ' 3 er g  1 ' 2 ) ,  and  the  fifth  number  is  the 
polarizability  (in  cubic  Angstroms).  These  data  were  taken  from  Kee,  et 
al.,13  except  for  the  potassium-containing  molecules.  The  transport 
parameters  for  K,  KOH  and  KO  were  assumed  to  be  like  those  of 'their  sodium 
analogs.  The  transport  parameters  for  KO^  were  assumed  to  be  like,  those  of 
Oj,  while  those  of  KH  were  assumed  to  be  like  argon.  For  each  species  there 
then  follow  the  14  thermodynamic  parameters  for  use  in  the  polynomial  fits  of 
Gordon  and  McBride.  These  parameters  are  tabulated  for  all  species  we 
consider  except  K02.  The  derivation  of  the  K02  parameters  is  discussed  in  the 
appendix.  The  use  of  these  parameters  is  detailed  elsewhere.  4 

Next,  we  limited  our  study  to  the  H2/air  flame.  The  choice  of 
fuel/oxidizer  ratios  was  somewhat  arbitrary  except  that  we  wanted  to  encompass 
a  reasonably  bread  range.  Figure  1,  taken  from  reference  17,  shows  the  flame 
velocity  of  H2/air  mixtures  vs.  H2  mole  fraction.  A  mole  fraction  of  0.25  for 
H2  gives  a  lean  mixture  while  one  of  0.50  gives  a  rich  one.  The  flame 
velocities  as  read  from  this  figure  are  182.5  cm/sec  and  287.5  cm/sec, 
respectively.  In  general,  a  given  value  of  the  flame  velocity  corresponds  to 
two  values  of  the  mole  fraction  and  we  simply  took  the  other  two  H2  mole 
fractions  for  our  test  cases.  This  procedure  gave  us  values  of  the 
equivalence  ratio  ranging  from  0.79  to  3.88,  as  can  be  seen  in  Table  7. 

This  table  shows  the  mole  fraction  initial  conditions  for  the  fuel,  H2, 
the  oxidizer,  02,  the  diluent,  N2,  and  the  inhibitor/suppressant,  expressed  as 
KOH.  It  also  sliows  the  computed  flame  speeds,  S(cm/s),  and  the  adiabatic 
flame  temperature,  T(Ke!vin).  : 

13.  J.  M.  Heimer  1  and  T.  P.  Coffee.  "The  Detailed  Modeling  of  Premixed, 
Laminar  Steady-State  Flames.  I.  Ozone,”  Combustion  and  Flame,  Vol.39,  pp.301- 
315,  1980. 

14.  T.  P.  Coffee  and  J.  M.  Heimerl,  "Transport  Algorithms  for  Premixed, 
Laminar  Steady-State  Flames,"  Combustion  and  Flame.  Vol.43,  pp.273-289,  1981. 

15.  R.  J.  Kee,  J.  Warnatz  and  J.  A.  Miller,  "A  FORTRAH  Computer  Code  Package 
for  the  Evaluation  of  Gas-Phase  Viscosities,  Conductivties,  and  Diffusion 
Coefficients,"  Sandia  Report  SAND83-8209,  March  1983. 

16.  S.  Gordon  and  B.  J.  McBride,  "Computer  Program  for  Calculation  of  Complex 
Chemical  Equilibrium  Compositions,  Rocket  Performance,  Incident  and  Reflected 
Shocks,  and  Chapman-Jouguet  Detonations,"  NASA-SP-273,  1971.  (1982  program 
version. ) 

17.  J.  Warnatz,  "Calculation  of  the  Structure  of  Laminar  Flat  Flames  II: 
Flame  Velocity  and  Structure  of  Freely  Propagating  Hydrogen-Oxygon  and 
Hydrogen-Air  Flames. "Ber.  Bunsenges.  Phvs.  Chem..  Vol.82,  pp.643-649,  1  17S. 
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TABLE  6.  TRANSPORT  AND  THERMODYNAMIC  PARAMETERS 


H  1.  2.05  145.  0.  0. 

.25000000D+01  0.  0.  0.  0.  H 

.25474390D+05  -.45989841D+00  .25000000D+01  0.  0.  H 

0.  0.  .254743 90 D+05  -.459898410+00  H 


OH  17.  2.75  80.  0.  0. 

.28897814D+01  .10005879D-02  -.22048807D-06  .201912880-10  -.39409831D-15  OH 

.38857042D+04  .555664270+01  .38737300D+01  -.13393772D-02  .16348351D-05  OH 

-.52133639D-09  .41826974D-13  .35802348D+04  .342024060+00  '  OH 


0  16.  2.75  80.  0.  0. 

. 2534296 1D+01  -.12478170D-04  -.12562724D-07  .690298620-11  -.637970950-15  0 

.292311080+05  .496285910+01  .303094010+01  -.22525853D-02  .39824540D-05  0 

-.32604921D-08  .10152035D-11  .291365260+05  .26099342D+01  0 

H02  33.  3.46  107.  0.  1.6 

.401730600+01  .22172883D-02  -.577101710-06  .713728820-10  -.36458591D-14  H02 

-.114124450+04  .378460510+01  .359641020+01  .52500748D-03  .75118344D-05  H02 

-.95674952D-08  .36597628D-11  -.893335020+03  .66372671D+01  H02 

H2  2.  2.92  38.  0.  .79 

.30558123 D+01  .59740400D-03  -.16747471D-08  -.21247544D-10  .25195487D-14  H2 

-.861684760+03  -.172070730+01  .294323270+01  .34815509D-02  -.77713819D-05  H2 

•74997496D-08  -.25203379D-11  -.976954130+03  -.18186137D+01  H2 

02  32.  3.372  128.7  0.  1.6 

.361221391*01  . 74853 166D-03  -.19820647D-06  .33749008D-10  -.23907374D-14  02 

-.11978151D+04  .36703307D+01  .37837135D+01  -.30233634D-02  .99492751D-05  02 

-.98189101D-08  .33031825D-11  -.10638107D+04  .36416345D+01  02 


H20  18.  2.6  572.  1.844  0. 

.27167633D+01  .29451374D-02  -.80224374D-06  .102266820-09  -.484721450-14  HjO 
-.29905826D+05  .663056710+01  .407012750+01  -.110844990-02  .415211800-05  H20 
-.29637404D-08  .807021030-12  -.302797220+05  -.322700460+00  HjO 


N2  28.  3.62  97. r  0.  1.76 

.28532899D+01  .160221280-02  -.629368930-06  .114410220-09  -.780574650-14  N2 

-.890080930+03  .639648970+01  .370441770+01  -.142187530-02  .286703920-05  N2 

-.12028885D-08  - . 1 39546 77D-1 3  -.10640795D+04  .223362850+01  N2 


K02  71.  4.69  39.4  0.  0. 

0.56193353D+01  0. 1550941 7D-02-0.68222055D-06  0.13299869O-09-0.95542863O-14 
-0.14355891D+05-0. 17198447D+01  0.34337129D+01  0.67246085D-02-0.33216077D-05 
-0.16047227D-08  0.13795100D-11-0. 137206890+05  0.97559016D+01  FIGGER, ANDREWS 

K  39.  3.567  1375.  0.  0.  LIKE  Na 

0.25673650D+01-0. 149335960-03  0.12342444D-06-0.53394240D-10  0.11948426D-13 
0.99550531D+04  0.466420810+01  0.2493C967D+01  0.501641770-04-0. 127512240-06 
0. 135404910-09-0. 51145936D-13  0.99786360D+04  0.50560438D+01  TRANDATA2 
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TABLE  6.  TRANSPORT  AND  THERMODYNAMIC  PARAFETERS( continued) 


KOH  56.  3.804  1962.  0.  0.  LIKE  NaOH 

0.56400949D+01  0.12510226D-02-0.34984547D-06  0.44566993D-1 0-0. 20870279D-14 

-0.29698732D+05-0.40 568187D+01  0.40733441D+01  0. 97217945B-02-0. 15988804D-04 
0. 1 21483 53D-07 -0.33 709342D-1 1-0. 2 9506 55 8D+05  0.292223738*01  TRANDATA2 

KO  55.  3.812  383.  0.  0.  LIKE  NaO 

0.44244778D+01  0.19936155D-03-0.37128837D-07  0.71308300D-11-0.50369687D-15 
0. 7205233 1D+04  0.32945055D+01  0.37410778D+01  0. 312420170-02-0. 48020039D-05 
0.34660605D-08-0.93599791D-12  0. 733687 14D+04  0.655376090*01  TRANDATA2 

KH  40.  3.542  93.3  0.  0.  LIKE  Ar 

0.39603386D+01  0.72190323D-03-0.26918715D-06  0.526173000-10-0.37872683D-14 
0.13501837D+05  0.84218060D+00  0.28157756D+01  0. 398710600-02-0. 33410548D-05 
0.88602942D-09  0. 11402847D-12  0.13805838D+05  0.671201450*01  TRANDATA2 


TABLE  7.  INITIAL  CONDITIONS  FOR  AND  SELECTED  OUTPUT  FROM 
THE  STEADY-STATE,  LAMINAR,  FLAT-FLAME  MODEL 


INITIAL  CONDITIONS  OPTPUT 

equivalence 


rat.!? 

XH2 

X02 

xn2 

XKOH 

1 

I 

0.79 

.25000 

.15750 

.59250 

.0 

174.6 

2168 

.24950 

.15740 

.59210 

.001 

182.3 

2164 

.24875 

.15671 

.58954 

.005 

152.1 

2153 

1.28 

.35000 

.13650 

.51350 

.0 

285.4 

2338 

.34965 

.13636 

.51299 

.001 

315.6 

2334 

.34825 

.13582 

.51093 

.005 

335.7 

23  20 

.34650 

.13514 

.50836 

.01 

296.9 

2302 

.34300 

.13377 

.50323 

.02 

192.6 

2267 

2.38 

.50000 

.10500 

.39500 

.0 

275.1 

1938 

.50000 

.10500 

.39426 

.00074 

297.9 

1937 

.49750 

.10448 

.39302 

.005 

280.2 

1920 

.49500 

.10395 

.39105 

.01 

132.6 

1902 

3.88 

.62000 

.07980 

.30020 

.0 

177.0 

1581 

.61938 

.07972 

.29990 

.001 

177.6 

1577 

.61876 

.07964 

.29960 

.002 

131 .0 

1574 

An  analysis  (a  so-called  screening  analysis)  of  each  of  these  flames 
shoved  that  several  other  reactions  could  be  eliminated.  These  reactions  were 
found  to  contribute  lejs  than  1Z  to  the  formation  or  destruction  of  any 
species  over  the  computational  space.  The  results  of  these  eliminations  is 
the  32-reaction  network  which  is  shown  in  Table  8.  (The  reactions  retain  the 
numbering  of  Table  5.) 


>' 
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TABLE  8. 

ABBREVIATED  REACTION  NETWORK 

7  7 

reaction 

A 

I 

c 

/ 

1  OH 

+  h2 

-  h2o 

+  H 

0.17E-15 

-1.6 

-3298.8 

/  ■ 

2  H 

♦  02 

*  OH 

♦  0 

0.20E-06 

0.9 

-16573.2 

3  0 

♦  H2 

-  OH 

♦  H 

0.25E-16 

-2.0 

-7551.4 

4  H 

♦  02  +  M 

-  ho2 

♦  M 

0.55E-29 

0.8 

0.0 

5  H 

+  ho2 

-  OH 

+  OH 

0.25E-09 

0.0 

-1003.5 

6  H 

♦  HO2 

-  h2 

♦  o2 

0.42E-10 

0.0 

-693.5 

7  OH 

+  H02 

-  h2o 

+  o2 

0.33E-10 

0.0 

0.0 

1 

8  0 

+  ho2 

-  OH 

♦  0, 

0.33E-10 

0.0 

0.0 

9  H 

+  H  +  M 

-  h2 

♦  H 

0.50E-29 

1.0 

0.0 

10  H 

♦  H  ♦  H2 

-  H2 

♦  h2 

0.27E-30 

0.6 

0.0 

11  OH 

+  OH 

-  0 

+  h2o 

0.25E-14 

-1.1 

0.0 

12  H 

K02 

-  KO 

+  OH 

0.12E-09 

0.0 

0.0 

13  KO 

♦  h2° 

-  OH 

+  KOH 

0.17E-09 

0.0 

0.0 

1 

14  H 

KOH 

-  K 

♦  h2o 

0.17E-10 

0.0 

0.0 

16  H 

*  ko2 

-  0 

♦  KOH 

0.33E-11 

0.0 

0.0 

17  H 

♦  ko2 

-  K 

♦  ho2 

0.33E-10 

0.0 

0.0 

19  K 

♦  ho2 

-  KO 

♦  OH 

0.33E-11 

0.0 

0.0 

21  0 

♦  ko2 

-  KO 

♦  o2 

0.17E-11 

0.0 

0.0 

22  H, 

♦  KO 

-  H 

♦  KOH 

0.83E-09 

0.0 

0.0 

| 

/ 

28  K 

+  02  ♦  M 

-  ko2 

♦  M 

0.83E-30 

0.0 

0.0 

29  KO 

♦  KO 

•  K 

♦  ko2 

0.17E-09 

0.0 

0.0 

33  H 

♦  KH 

-  K 

+  h2 

0.17E-09 

0.0 

0.0 

38  KH 

♦  OH 

-  K 

♦  H20 

0.17E-09 

0.0 

0.0 

39  KH 

♦  O2 

-  K 

♦  ho2 

0.17E-09 

0.0 

0.0 

1 

40  H 

♦  ko2 

»  KH 

♦  o2 

0.17E-09 

0.0 

0.0 

42  CO 

+  OH 

-  H 

♦  co2 

0.73E-17 

-1.5 

741.0 

43  CO 

♦  o2 

♦  HO, 

»  0 

♦  C02 

0.42E-11 

0.0 

-47800.0 

44  CO 

-  OH 

co2 

0.25E-09 

0.0 

-23590.0 

45  CO, 

♦  0 

-  CO 

+  Oj 

0.28E-10 

0.0 

-52580.0 

46  CO 

♦  H 

.  CO 

♦  OH 

0.27E-09 

0.0 

-26300.0 

47  0 

♦  CO  +  M 

-  CO, 

♦  M 

0.15E-33 

0.0 

4540.0 

i 

48  KO 

♦  CO 

«  K 

♦  co2 

0.17E-09 

0.0 

0.0 

1 

The  rate  coefficient  ■  A  x 

r»  «  .C/,T 

where  A  in 

units  of 

cm-particle-sec 

C  in 

units  of 

calorie/mole 

i 

This  32-reaction  network  was  also  used 

as  input  for  the  MEFF  code 

computations  and  results  similar  to 

the  48  reaction  network  were  obtained.  In 

addition, both  hydrogen-air  subsets  of  Tables  5 

and  8, 

when  used  as  input  to 

the  flat-flame  code,  yielded  nearly 

identical  species. 

temperature,  and  heat 

1 

release 

profiles. 

We  conclude  that 

the  abbreviated  reaction  network  of  Table 

8  constitutes  a  minimal  descriptive 

suppression  network. 

1 

Caution  must 

be  urged 

at  this 

point.  If  the 

values 

of  any  of  the  assumed 

rate  coefficients, 

upon  measurement, 

prove  to  be 

significantly  less  than  given 

in  Table  8,  then,  because  of 

the  non 

-linear  couplings  among  all  of  the  species 

and  between  the  species  and 

the  temperature,  the 

entire 

network  given  in  Table 

1 

5  should  be  reexamined. 

! 

i 

/ 
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Referring  to  Table  7,  one  can  see  that  the  computed  flame  velocity  drops 
with  continued  addition  of  suppressant,  while  the  temperature  is  hardly 
affected.  This  is  not  surprising  since  we  are  dealing  with  a  series  of 
different  flames  and  there  need  not  be  any  simple  correlation  among  them.  To 
see  this  most  easily,  consider  the  two  flames  with  initial  fuel  and  oxidizer 
values  of  0.50000  and  0.10500,  respectively  (equivalence  ratio  *  2.38).  The 
only  difference  lies  in  the  trace  amount  of  KOH  added  at  the  expense  of  the 
diluent  No*  The  adiabatic  flame  temperature  changes  by  one  degree  but  the 
computed  flame  speed  shows  a  noticeable  change.  Analysis  shows  that  near  the 
peaks  of  the  respective  heat  release  functions,  the  elementary  chemistry  has 
altered  radically.  In  the  case  of  no  KOH  addition,  the  rate  of  increase  of 
temperature  is  governed  by  the  reactions: 

H  +  H02  -  OH  +  OH,  and 

H  +  H02  -  H2  +  02, 

and  by  diffusion  processes.  Iu  the  case  of  the  small  addition  of  KOH,  the 
rate  of  increase  of  temperature  is  governed  by: 

H  +  02  +  M  -  H02  ♦  M, 

K  ♦  02  ♦  M  -  K02  +  M, 

OH  +  H2  ■  H20  +  H,  and 

H  H02  *  OH  ♦  OH; 

all  listed  in  rank  order. 

The  above  points  out  the  importance  of  obtaining  reliable  values  for  key 
reactions  such  as  K  +  02  +  M  *  K02  ♦  M.  Were  we  to  repeat  any  of  these 
computations  we  would  use  the  most  recent  measurements  of  Husain  and  Plane18 
and  of  Silver  et  al.  ®  Their  values  for  this  potassium  reaction  show  that  our 
assumed  value  is  low  by  a  factor  ranging  from  5  at  300  K  to  5QZ  at  3000  K. 
Had  we  used  this  value  somewhat  more,K02  would  have  been  formed;  however,  this 
should  not  qualitatively  affect  our  results. 

The  fact  that  the  rate  coefficient  for  the  reaction  forming  K02  is  so 
great  has  implications  for  the  experimentalist.  When  we  first  began  to  use 
the  flat-flame  code  we  initialized  all  of  the  potassium  as  the  atomic  species, 
K.  Analysis  showed  that  within  the  space  scale  we  were  using  (about  0.01  cm), 


18.  D.  Husain  and  J.  M.  C.  Plane,  "Kinetic  Investigation  of  the  Reaction 
between  K  +  02  +  M  by  Time-Resolved  Atomic  Resonance  Absorption  Spectroscopy," 
J.  Chem.  Soc.  Faraday  Trans.  2.  Vol.78,  pp.1175-1194,  1982. 

19.  J.  A.  Silver,  M.  S.  Zahniser,  A.  C.  Stanton,  and  C.  E.  Kolb,  "Temperature 
Dependent  Termolecular  Reaction  Rate  Constants  For  Potassium  and  Sodium 
Superoxide  Formation,"  Aerodyne  Report  AR1-RR381,  January  1984.  See  also: 
paper  #68  Twentieth  Symposium  (International)  on  Combustion,  the  University  of 
Michigan,  Ann  Arbor  Michigan,  12-17  Aug  1984. 
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all  Che  K  was  converted  to  K02  at  the  cold  (298  K)  boundary.  The  conclusion 
we  draw  is  that  unless  special  precautions  are  taken,  the  addition  of  free 
atonic  potassium  to  an  input  gas  stream  containing  02  w°uld  result  in  the 
conversion  of  the  K  to  K02  before  the  gases  could  reach  the  burner  surface. 
(In  the  execution  of  the  flat-flame  computer  code  this  difficulty  was  easily 
overcome  by  assigning  all  the  initial  potassium  to  the  species  KOH.) 

Table  8  shows  that  the  species  KH  is  involved  in  only  a  few  reactions. 
Analysis  shows  that  the  KH  is  almost  exclusively  formed  by  the  reaction: 

H  +  K02  -  KH  ♦  02. 

The  computed  concentration  of  K02  can  be  large  for  an  intermediate  species  and 
in  our  computed  results  it  is  probably  too  low.  To  be  able  to  tell  whether  or 
not  the  species  KH  will  actually  play  a  significant  role  in  the  flash 
suppression  network  will  require  a  measurement  of  this  reaction  rate 
coefficient . 


IV.  SUMMARY 


An  advanced  suppression  network  has  been  constructed  by  extending  the 
inhibition  network  of  Hynes  et  al.^  First,  the  most  likely  species  to  be 
involved  in  flash  suppression  were  considered.  These  species  include:  K,  KO, 
K02,  KOH, and  KH.  From  these  species  a  complete  set  of  55  two— body  and  15 
three-body  reactions  was  constructed.  Standard  chemical  kinetic  arguments 
were  used  to  reduce  the  number  of  reactions  to  48.  Use  of  screening  analysis 
upon  the  output  of  flat-flame  code  computations  enabled  the  network  to  be 
further  reduced  to  32  reactions.  Both  the  extended  set  of  48  reactions  and 
the  abbreviated  set  of  32  reactions  gave  nearly  identical  results  when  used  as 
the  reaction  network  for  the  MEFF  muzzle-flash-prediction  code,  for  both 
howitzer  and  mortar  calculations.  In  addition,  when  used  as  input  for  the 
flat-flame  code,  the  respective  hydrogen-air  subsets,  with  and  without 
potassium,  also  yielded  the  same  computed  profiles  for  Che  species, 
temperature,  and  heat  release.  Further  analysis  showed  that  K02  is  mainly 
formed  from  K  ♦  02  ♦  M  *  K02  ♦  M  and  that  KH  is  mainly  formed  from 
H  +  K02  »  KH  ♦  0^.  The  former  is  currently  under  investigation1,9  while  the 
latter  is  a  candidate  for  future  work. 
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APPENDIX  A 


Approximate  Properties  Of  The  KO2  Molecule 


The  one-dimensions  1 ,  laminar,  flat-flame  code  we  used  demands  species- 
dependent  thermodynamic  functions  of  temperature.  This  was  accomplish  by 
using  the  polynomial  expressions  of  Gordon  and  McBride.^-*  For  each  reaction 


species,  the  specific  heat,  enthalpy,  and  entropy  as  functions  of  temperature 
are  given  in  the  form: 


Cp/R  -  ai  +  a2T  +  a3T2  +  a^T3  +  , 

Ht/RT  -  ax  +  a2T/2  +  a3T2/3  +  a^M  +  a5T4/5  +  a6/T, 

and 

Sf/R  m  aj InT  +  a2T  +  a3T2/2  +  a^T3^  +  a^T^/4  +  a 


Parameters  corresponding  to  each  species  considered  in  this  report  are 
documented  except  for  the  species  K02.  Suitable  parameters  for  K02  can  be 
generated  provided  that  the  heat  of  formation,  moments  of  inertia,  and  the 
normal  mode  frequencies  are  known.  First  consider  an  estimation  of  the  K02 
heat  of  formation. 


Figger,  Schrepp,  and  Zhu^-2  used  spectroscopic  techniques  to  determine  a 
lower  limit  for  the  K-02  bond  enery  of  (45.2  +/-  2.4)kcal/mole.  It  is  our 
present  task  then  to  relate  this  bond  energy  to  the  heat  (enthalpy)  of 
formation. 


The  standard  enthalpy  of  formation  of  any  compound  is  the  net  heat  of 
reaction  by  which  it  is  formed  from  its  elements,  the  reactants  and  products 
all  being  in  a  given  standard  state.  A  convenient  standard  state  for  a 
substance  ij  the  state  in  which  it  is  stable  at  298. 15K  and  one  atmosphere 
pressure.  By  convention  the  enthalpies  of  the  chemical  elements  in  this 
particular  state  are  set  equal  to  zero. 

Using  the  first  law  of  thermodynamict  and  considering  only  PV  work,  we 
find  the  change  in  enthalpy  can  be  written  as: 

''  AI H  -  4.U  +  &.(PV). 

By  A_(PV)  we  mean  the  change  in  PV  for  the  entire  system;  that  is,  the  PV  of 
the  products  minus  the  PV  of  the  reactants  for  the  particular  chemical 
reaction.  For  reactions  in  the  gas  phase,  the  values  of  depend  on  the 

change  in  the  number  of  moles  of  gas  as  a  result  of  the  reaction  under 
consideration.  From  the  ideal  gas  equation  we  can  write: 

•L(PV)  =  RTAn. 


A-l.  S.  Gordon  and  B.  J.  McBride,  "Computer  Program  for  Calculation  of 
Complex  Chemical  Equilibrium  Compositions,  Rocket  Performance,  Incident  and 
Reflected  Shocks,  and  Chapman-Jouguet  Detonations,"  NASA-SP-'273,  1971.  (1982 
program  version.) 

A-2.  H.  Figger,  W.  Schrepp  and  X.  Zhu,  "Chemiluminescent  Reaction  Between 
Alkali  Dimers  and  Oxygen  Molecules."  J.  Chem.  Phvs..  Vol.79,  pp  1320-1325, 
1983  . 
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Thus : 


A_h  -  iU  +  RT~n, 

where  by  An  we  mean  the  number  of  moles  of  gaseous  products  minus  the  number 
of  moles  of  gaseous  reactants. 

For  the  present  case: 

K02  =  K  +  02- 

We  note  that  at  zero  degrees  absolute: 

AHC  -  AUq  -  ^ r- 

Now : 

£Hr  =  AHf0(K)  +  AH£o(02)  -  AHfo(K02) 
and  substituting  values  from  Table  A-l  we  find 
45.2  -  21.5  +  0.0  -  £Hfo(K02) 
or 

£_Hfo(K02)  ■  -23 ,7kcal/mole. 


TABLE  A-l.  SELECTED  THERMODYNAMIC  VALUES 


Species 

— f298 

±tz. 

(H(pH298l 

Referen 

K 

21.522 

21.31 

0.20 

-1.481 

A-8 

°2 

0.0 

0.0 

-2.075 

A- 8 

°3 

34.8 

34.2 

0.4 

-2.474 

A-8 

h2o 

-57.103 

-57.7979 

-2.367 

A-8 

KOH 

-54.6 

-55.6 

3.0 

-2.730 

* 

ko2 

(bond  energy) 

-  (45.2 

+/-  2 

.4)  kcal/mole 

A- 2 

*  J.  Chem.  Phvs.  Ref.  Data.  Vol.3,  p.432,  1974. 


To  find  the  heat  of  formation  at  298. 1  5K  we  recognize  that  for  the  reaction 
above  where  the  stoichiometric  coefficients  are  unity: 

-HfT  «■  ^_Hf298  +  (Ht  -  H29g)compound  -  SUM(Ht  -  H29g), 

where  the  sum  is  taken  over  the  elements  in  their  standard  state.  For  this 
case  we  have: 

4_Hf()(K02)=  AHf298(K02^  +  ^H0-H298V.02  '  ^H0“H298^K  “  ^Ho"H298^02* 
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Substituting  appropriate  values  from  Table  A-l  we  find: 

-Uf298^K02^  *  “(25  ♦/-  3)  kcal/mole, 

where  we  have  approximated  the  difference  in  enthalpy  for  K02  by  that  for  O3. 

As  a  check  on  this  value  we  note  that  the  room  temperature  heat  of 
reaction  for: 

303  +  2K0H  -  302  ♦  2K02  +  H20  -  101.2  kcal/mole, 

was  reported  by  the  Vol'nov  et  al  .4-3  Using  the  values  given  in  Table  A-l,  we 
find: 


^.Hf2^g(K02)  ■  -(26  +/-  3)  kcal/mole, 
in  good  agreement  with  our  above  estimate. 

Both  Raman  spectra^-^  and  electron  spin  resonance  spectraA”^  of  K02  in 
rare  gas  matrices  show  that  the  structure  is  that  of  an  isosceles  triangle. 
Since  the  actual  structure  of  the  gas  phase  molecule  has  not  been  determined, 
we  assume  that  the  gas  phase  molecule  also  has  an  isosceles  triangle 
structure,  identical  with  that  of  the  solid.  Andrews11”®  estimated  the  K-0 
distance  in  K02  to  be  2.28A,  while  Kassatochkin  and  Kotow^*"  obtained 
(1.28  +/-  0.07)A  as  the  distance  between  adjacent  oxygen  atoms.  With  these 
data  the  moments  of  inertia  can  be  computed.  They  are:  2.18,  14.0,  and  16.18 
in  units  of  10”-^gm-cm  . 


The  last  cats  we  require  for  the  K02  molecule  are  the  values  of  the 
normal  mode  fr«quencies.  Smardewski  and  AndrewsA”^  measured  the  fundamental 
frequency  to  be  1108  cm'  .  They  also  measured  the  symmetric  interionic 
stretching  mod  2  frequency  of  307.5  cm”*.  Based  on  analogy  with  ozoneA  we 
estimate  the  antisymmetric  normal  vibration  frequency  to  be  1100  cm”1. 
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Metal  Ozonide  •ormation."  Isv.  Akad.  Nauk.  5SSR.  Ser.  Khim,.  pp.1183-1187, 
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Noting  that  the  assumed  geometry  for  KO2  yields  a  symmetry  number  of  two,  we 
have  all  the  required  input  information  to  generate  the  Gordon  and  McBride 
coeff icients.  ’  The  values  for  the  a^  are  shown  in  Table  A-2. 


TABLE  A-2.  GORDON  AND  MCBRIDE  COEFFICIENTS,  ai(  FOR  R02 


T>1000K 


T<1000K 


1  .56193353E+01 

2  . 1550941 7 E-02 

3  -.68222055E-06 

4  .13299869E-09 

5  -.95542863E-14 

6  -.14355891E+05 

7  -.17198447E+01 


.34337129E+01 
.67246035 E-02 
-.33216077E-05 
-. 16047 227 E-08 
.13795100E-11 
-.13720689E+05 
•97559016E+01 


Use  of  these  a^  and  R*1.9872  cal/mole/K  yields  a  room  temperature  value 
for  the  entropy  of  K02  of  62  cal/mole/K. 


A-9.  B.  J.  McBride  and  S.  Gordon,  "FORTRAN  IV  Program  for  Calculation  of 
Thermodynamic  Data,"  NASA  IN  D-4097,  1967. 
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reports  it  publishes.  Your  comments/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ Date  of  Report _ - _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
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etc?  If  so,  please  elaborate. _ . _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
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